function P=IP_pos_ang(G,P,k)
P.pos=[P.Vars(1),0];
P.ang=P.Vars(3);

